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ABSTRACT 

Using wide-field, two-color HST imaging of the cluster CL1358H-62 (z = 0.33), we derive structural parameters 
for a large, magnitude-limited sample of cluster members. These structural parameters are combined with accurate 
velocity dispersions in another paper to investigate the fundamental plane in the cluster We fit integrated r'/"*- 
laws to the integrated surface brightness profiles, and fit two-dimensional r'/'^-law model galaxies to the images 
directly. A comparison of the results from the two methods shows that the derived half-light radii, r,,, agree very 
well, with an nns scatter of only 13%. The half-light radii range from approximately I to 20 kpc, with a median 
of about 3 kpc (Ho = 65kms"'Mpc"', = 0.1). 

We investigated the stability of r^, by comparing the r'/'^-law fits to the half-light radii derived using other 
commonly used surface brightness profiles. In particular, we fit Sersic r'/"-laws (finding the best-fit n ranging 
between n = 1-6) and superpositions of r'/'^-law bulges with exponential disks. The half-light radii derived from 
the best-fit Sersic profiles varied with respect to from the r'/'^-law fits by only 1% in the median, but with a 
standard deviation of more than 40% in . For the bulge-plus-disk fits, the derived half-light radii were offset 
from of the r'/'*-law fits by 10% in the mean, also showing a large standard deviation of more than 40%. By 
comparing the fitted half-light radii from the Sersic laws with those derived from the bulge-plus-disk fitting, one 
also finds a large scatter of 30%. Based on these tests, we conclude that, in general, half-light radii are only 
measured with a typical accuracy of 30-40%. 

While there are large uncertainties in half-light radii, these do not impact the subsequent fundamental plane 
analysis because the combination r{I)^-''^, which enters the fundamental plane, is extremely stable. The nns 
scatter in r{I)^-''^ is less than 3%, regardless of the form of the profile fit to the galaxies. 

By fitting bulge-plus-disk profiles, we find that the median bulge fraction of the sample is 84% and that the 
few Eh-A galaxies in this sample have disks which make up about 10-35% of their total light. These results are 
consistent with the residuals from fitting two-dimensional r'/'^-law models directly to the galaxies, in which disk- 
Uke structures are revealed in about half of the sample. Two of the three Eh-A galaxies show spiral arm structure. 



L INTRODUCTION 

Galaxy scaling relations provide a unique opportunity to 
study the evolution of stellar populatio ns through intermediate 
redshift s (iFranx 1993L IVogt ef a/. 19961 1997). Sinc e their dis- 
covery ( [Faber & Jackson 1976 , Fully & Fisher 1977 ), such scal- 
ing laws have been used primarily to measure distances, with a 
reliability seemingly limited only by measurement errors and 
the homogeneity of the galaxies in question. The Hubble Space 
Telescope (HST) and WFPC2 have enabled high S/N imag- 
ing of intermediate redshift galaxies with a resolution similar 
to ground-based imaging of galaxies in the Coma cluster Us- 
ing the high spatial resolution and S/N of WFPC2, and effi- 
cient spectrographs at large aperture telescopes like the Keck 
10m, we can directly measure galaxy scaling relations up to 
redshifts of unity in an ef fort to understand t heir origins, and 
how gala xies evolve (e.g., ISchade et al. 1997 , van Dokkum et 
al. 1998b, |Simardefa/.1999D . 

Of the k nown scaling relations, the fundamental plane (F aber 
et al. 1987, pjorgovski & Davis 1987[ ) has the lowest observed 
intrinsic scatter, only -15% in dist ance for a given galaxv (J 0r- 
gensen, Franx, & Kjaergaard 1993, Kelson et al. 1997). The 



fundamental plane (FP) is an empirical correlation between ef- 
fective radius, r^, central velocity dispersion, a, and surface 
brightness, , for early-type galaxies. In a study of 1 1 nearby 
clusters J0rgensen et al. (1996) found, in Gunn g. 
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with a low (rms) scatter of ±20% dex in . In Coma, the scat- 
ter i s even lower, ±14% in distance (Luce y, Bower, <& Ellis 
1991, [l0rgensen, Franx, & Kjaergaard 1993 ). The fundamental 
plane of Virgo elli pticals also has a v ery low internal scatter of 
±12% in distance ( Kelson era/. 1999 ). Using the virial theorem 
and assuming that early-type galaxies are an homologous fam- 
ily, the low scatter in the fundamental plane implies a scatter of 
only ±20% in M/L ratio for a given a galaxy mass. 

The utility of the fundamental plane for probing the lumi- 
nosity evolution of early-type galaxies has been clearly demon- 
strated by van Dokkum & Franx (1996), Kelson et al. (1997), 
and van Dokkum et al. (1998b), who measured the evolution of 
the fundamental plane zero-point out to z = 0.83. The implied 
evolution of the MjLy ratios is consistent with a high redshift of 
formation for stars in early-type galaxies in rich clusters. Pro- 
jections of the fundamental plane, such as the r-L relation (Ko- 
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rmendy 1977) and the 1 -parameter L-a Faber- Jackson relation 
(Faber & Jackson 1976), have also been used by other authors 
to study e arly-type galaxy evolution to high redshifts TS chade 
et al. 1997, |Zie.gler & Bender 1997L [Bender et al. 1998|, Pahre, 
Djorgovski, & Carvalho 1998). 

Other studies of high redshift cluster galaxies have confirmed 
the notion t hat early-type galaxies are homogeneous , and likely 
to be old dSchade et a/. 19971 lElHs et al. 19971). Further- 



more, Dressier et al. (1997) showed that there has been substan- 
tial evolution in the numbers of SO galaxies between redshifts 
z ^ 0.5 and the present. However, many such results derive 
from limited imaging of cluster cores. Recently, van Dokkum 
et al. (1998a) used the color-magnitude relation of early-type 
galaxies to show that the uniformity of age extends to large radii 
for the ellipticals of CL1358H-62, and that the most recent star- 
formation in massive clusters has occurred in the outer parts, in 
galaxies with disks, van Dokkum et al. showed that the stellar 
populations of galaxies with disks are less homogeneous than 
pure bulge systems, but only at cluster radii which are outside 
the standard WFPC2 footprint (at the distances of intermedi- 
ate redshift clusters). Evolution in cluster galaxies may be best 
detected in large, wide-field surveys of clusters: cluster mem- 
bership, wide spatial extent, and large samples inclusive of all 
morphological types are all crucial for finding, and measuring 
the evolution in cluster populations. Surveys which are specif- 
ically biased towards ellipticals may miss much of the recent 
(z < 1) evolution in cluster galaxies because the samples are 
biased against those morphologies which have undergone the 
most recent epochs of star-formation. 

The first detailed study of membership in CL1358H-62 was 
performed by Fabricant et al. (1991), who measured a blue 
and active galaxy fraction of -15%. Complete membership to 
R = 21 mag is discussed by Fisher etal. (1998) and a small sub- 
sample was used by Kelson etal. (1 997) to find the fundamental 
plane zero-point for the cluster van Dokkum et al. (1998a) has 
made a complete analysis of the Color-Magnitude relation of 
194 confirmed cluster members within the 8' x 8' HST mosaic. 
The imaging consists of two colors (F606W and F814W) in 12 
CVZ pointings. This WFPC2 imaging of CL1358H-62 is a su- 
perb foundation for building a detailed history of cluster galax- 
ies, from the current epoch, through intermediate redshifts. Fu- 
ture work will concentrate on mosaics of WFPC2 imaging in 
the clusters MS2053-04 (z = 0.58) and MS 1054-03 (z = 0.83). 

Since the publication of Kelson et al. (1997), we have com- 
pleted the surface photometry for a high-resolution spectro- 
scopic sample of galaxies within the HST mosaic's field. This 
subset of the large membership database will be used for analy- 
sis of the fundamental plane and absorption line strength prop- 
erties. The data are discussed in §^. The photometric cali- 
bration, the transformation to redshifted Johnson colors, and 
sources of systematic calibration errors, are discussed §|^. The 
techniques used to fit r'''''-law profiles to the data are discussed 
at length in §^. In §|| we discuss such instrumental issues as the 
under-sampling and point-spread function errors. Deviations 
from the de Vaucouleurs profile a re explored a t some length in 
Sections ^ and fusing r'/"-laws (Sersic 1968) and bulge-disk 
decomposition of the circularized profiles. Tables of the struc- 
turalparameters are given in §^and our results are summarized 
in §1 

2. THE DATA 

The HST imaging consists of 12 contiguous pointings, mak- 
ing up a large mosaic of the cluster environment, about 1.5/i~' 



Mpc on a side (where Hq = lOO/z kms~'Mpc~'). Each point- 
ing consists of three 1200 s exposures in each of F606W and 
F814W. The data were preprocessed by the standard pipeline 
proc edures which are discussed elsewhere ( Holtzman et al. 
1995). 

The analysis presented in this paper was performed on the 
images processed by van Dokkum et al. (1998a), who cleaned 
the images of cosmic rays, modeled and removed the compli- 
cated background, and combined the individual exposures. The 
most demanding task was the removal of the scattered earth 
light that resulted from the observations being performed in 
the C ontinuous Viewing Zone (CVZ) (see v an Dokkum et al. 
1998a). In some of these observations, the background was 
higher by factors of two to three. This scattered light increases 
the counts in the sky background, and the secondary mirror 
supports cast shadows in the pattern of an "X," in which the 
background counts are typically not elevated. The excess in 
the background, and contrast between the shadowed and illumi- 
nated regions of the CCDs can vary greatly, and its magnitude 
depends on the position of the telescope in its orbit when a CVZ 
exposure is taken. Fortunately, in each set of three exposures, 
in a given pointing and filter, at least one frame had minimal 
or no scattered light, thereby allowing accurate modeling of the 
scattered earth light. Large-order (n i 15) bivariate polynomials 
were fit to the background and subtracted. More details of this 
modeling are given by van Dokkum et al. (1998a). 

The selection of our fundamental plane sample is discussed 
in Kelson et al. (1999a). Briefly, we randomly selected clus- 
ter members within the field of view of the HST mosaic, to 
R <21 mag. The selection was performed to efficiently con- 
struct mu hi-slit plates for the Low-Resolution Imaging Spec- 
trograph (Oke et al. 1995). We used three masks, with different 
position angles on the sky. The region of maximum overlap is 
in the center of the cluster, and thus the FP sample is concen- 
trated towards the core of the cluster 

Morphology was not a factor in the selection of our sample. 
In the random selection process, three E+A galaxies were in- 
cluded; these will be compared with those cluster members that 
have normal, early-type spe ctra. The sample is about 50% com- 
plete for R < 20.5 mag (see Fisher et al. 1998 for details on the 
statistical completeness of the original redshift catalog). Three 
galaxies fainter than R = 2l mag were added to test the quality 
of velocity dispersion measurements for faint galaxies. These 
spectroscopic observations and their reductions are detailed in 
Kelson et al. (1999a). In total, we have central velocity dis- 
persions for 55 galaxies in the cluster Unfortunately, two of 
the galaxies were imaged too close to the WFPC2 CCD edges 
to obtain reliable structural parameters. Thus, there are a total 
of 53 cluster members in our determination of the fundamental 
plane in CL1358H-62. 

Images of these 53 galaxies are shown in Figure |l} The im- 
ages clearly show a wide range of morphologies, from ellip- 
ticals to early- type spirals. Optical morphologies have been 
taken from Fabricant et al. (1999). These authors classified 
several hundred galaxies in the HST mosaic according to mor- 
phological type, T. The galaxies in our high-resolution spec- 
troscopic sample have T e {-5,-4,-3,0, 1,2,3} (E, E/SO, SO, 
so/a, Sa, Sab, Sb, respectively). The galaxy morphologies are 
listed in Table ^ in which there are 11 E, 6 E/SO, 14 SO, 13 
so/a, 6 Sa, 2 Sab galaxies and 1 Sb galaxy. This distributio n is 
quite similar to nearby massive clusters (e.g., Oemler 1974). 



3. PHOTOMETRIC CALIBRATION 
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The aim of this work is to produce surface photometry for 
use in the analysis of the fundamental plane in the cluster To 
enable comparison of the CL1358+62 fundamental plane with 
equivalent observations of nearby galaxies, the calibrated HST 
photometry must be transformed to a standard system. Such 
transformations to a photometric system of redshifted John- 
son B and V bandpasses will allow us to compare our obser- 
vations directly to observations of nearby galaxies. Precisely 
these transformations and calibration steps have been used pre- 
viously by Kelson et al. (1997) and van Dokkumefa/. (1998a). 
Each step has its own uncertainties, which will be summarized 
at the end of the section. 

3.1. Zero-points 

In this section, we describe the multiple components in the 
zero-point calibration of the HST photometry in F606W and 
F8 14W. The starting points of the calibration are the zero-points 
for F606W and F8 14W derived by Holtzman et al. (1995) using 
photometry of stars (point sources) with an aperture of r = 0''5. 
Our measured surface brightnesses have an effective aperture 
defined by the extent of our adopted point-spread functions, in 
this case r = 1".5. Thus, we must correct our magnitudes from 
the r = l".5 aperture to an effective aperture of r = 0''5. Us- 
ing the encircled energy curves in the HST Data Handbook, we 
estimate the correction to be +0.067 mag. The Holtzman et 
al. (1995) calibration is also valid for data taken using a gain of 
14e"/ADU. Therefore we also incorporate the gain ratios for the 
individual CCDs as given by the handbook. Furthermore, the 
HST Distance Scale Key Project has measured a difference of 
0.050 ma g between long and short exposures of identical fields 
(see, e.g.Milletal. 19971, [Stetson ef aZ. 19981) such that long in- 



tegrations have more counts than would be expected by simply 
scaling the shorter exposures. Thus, the WFPC2 zero-points, 
which were derived from short integrations of star fields, are 
not appropriate for long integrations, and need to be corrected 
by +0.050 mag. 

In summary, the F606W and F814W magnitudes are defined 

by 

F606W = -2.5 log(DN/s) + 2 1 .419 + 2.5 log(GRi) + 0.067 + 0.05(Q) 
F8 14W = -2.5 log(DN/s) + 20.840 + 2.5 log(G/?,) + 0.067 + 0.0a;0) 

where GR, represents the ratios of the high- to low- 
gain modes for the four WFPC2 CCDs, respectively: 
1.987,2.003,2.006,1.955. The total exposure time for each 
point was 3600 s per filter 

3.2. Transformation to redshifted V, (B — V) 

In order to compare the surface brightnesses and M/L ratios 
of distant galaxies to equivalent observations made nearby, we 
must transform the photometry to a consistent system. Towards 
this end, we define B, and V, as the Johnson B and V filter band- 
passes redshifted to z = 0.33, the observed frame of the galax- 
ies. As one can see in Figure ^, the B, and V, bands fall at 
roughly 5500 A and 7000 A, overlapping substantially with the 
observed F606W and F814W bandpasses. Without the color 
information, one would be required to extrapolate to, e.g., V 
magnitude using K corrections. However, we can derive an ac- 
curate transformation from the WFPC2 system to the redshifted 
Johnson system by interpolation using the color information. 

As in van Dokkum & Franx (1996), the flux in the redshifted 
filter can be written as a function of the fluxes in the observed 



bandpasses: 

P'v, = -Ppeoew ^ -fraMW' (4) 
where a is a color term in this transformation (see below). In 
order to relate the flux in a given bandpass to a standard magni- 
tude, one must use filter constants (see below). Solving for as 
a function of the observed magnitudes and the filter constants, 
cx, finds: 



V- = F8 14W- CF8I4W + cy + fl [(F606W - F8 14W)- 
(cf606w-cf814w)] +2.51og(l+z) 



(5) 



The additional factor of 2.51og(l+z) compensates for the 
stretching in wavelength of the redshifted bandpass. The B^ 
and magnitudes are redshifted measures of the fluxes, not 
restframe quantities. Removal of the (l+z)'^ surface brightness 
dimming is required to convert these data to magnitudes in the 
restframe of the galaxies. 

The filter constants, cx, are used to convert the V,, F814W, 
and F555W magnitudes to a system of flux units in a given 
filter This system can be thought of as similar to an AB sys- 
tem. These constants are defined by cx = X-Xab, where X is 
the magnitude, through filter X, in the standard system. One 
can calculate the AB magnitudes of Vega in the different fil- 
ters, a nd, using the reported magnitudes in the literature for 
Vega (Holtzman et al. 1995), define the transformation from 



the standard (Johnson) system to an AB system, based on abso- 
lute flux uni ts. We used the Vega spectrum from Bohlin (1996, 
Hayes 198^ , the filter transmission curves for the WFPC2 fil- 
ters, and the Johnson filter curves (Bessell 1991), as well as the 



WFPC2 CCD response curve to determine the AB magnitudes 
of Vega in the various filters. These AB magnitudes are zero- 
pointed to Vega^ not the spectra of F dwarfs (i.e. AB79, O ke & 
Gunn 1983; see |Frei & Uunn 1994 , or Fukugita, Shimasaku, & 
Ichikawa 1995, Fukugita ef a/. 1996[ ). 

Using the SEDs of E/SO galaxies ( ^ence 1976 ), one inte- 
grates over the observed and redshifted filter curves to solve 
for a, for a given SED: 



(V,-cv)-(F814W-CF8i4w) 
(F606W - CF606W) - (F8 14W- cfsmw) 



(6) 



to find 



y_ = F814W+ 0.204 X (F606W-F814W) + 0.652 (7) 
(B,-y_) = 0.818 X (F606W-F814W)- 0.128 (8) 

for galaxies with E/SO spectra at z = 0.33. Using SEDs typi- 
cal of spirals yield transformations which differ by only ±0.01 
mag. As noted above, these transformations are identical to 
those used by Kelson et al. (1997) and van Dokkum et al. 
(1998a). 

3.3. Colors 

Colors in the WFPC2 filters are taken directly from van 
Dokkumefa/. (1998a). They measured (F606W-F814W) col- 
ors within for each of the cluster members within the HST 
mosaic; their color errors (random) are on the order of ±0.01 
mag. For the purposes of the photometric calibration, the abso- 
lute color uncertainties should approximately equal the quadra- 
ture sum of the two Holtzman et al. (1995) zero-point errors. 
Other effects, such as the exposure time correction, the charge- 
transfer effect, gain variations and aperture corrections, cancel 
since they are equal contributors to the calibration of each band- 
pass. We comment more on color uncertainties below. 
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3.4. Galactic extinction 

Galactic extinction has been taken from Burstein & Heiles 
(1982) and used in conjunction with Cardelli, Clayton, & 
Mathis (1989) to derive the extinctions in the HST filters. The 
Burstein & Heiles E(B-V) color excess is 0.006 mag. Us- 
ing Cardelli et ai, we find ^f606w = 2.91 and /?f8I4w = 1-94. 
Thus, Af6()6w = 0.02 and Apsww = 0.01. Thus, we correct both 
the measured colors and the redshifted V surface brightnesses 
by 0.01 mag. Because the Galactic reddening is low in the di- 
rection of CL1358H-62, reddening uncertainties will contribute 
little to the following error budget. 

3.5. Uncertainties in the Redshifted Photometry 

Errors in the final calibration of the redshifted y,-band sur- 
face brightnesses can arise from three sources: (1) errors in 
the HST zero-point calibration, (2) errors in the measured col- 
ors, and (3) errors in the transformation from the observed to 
redshifted bandpasses. We now attempt to summarize these 
uncertainties, and estimate the total error in our calibrated V~ 
magnitudes. 

Each step in the calibration of HST photometry introduces 
uncertainties. In the previous section, we showed that the V^- 
bandpass is very closely tied to F814W, with a moderate color 
term. Thus, errors in the Holtzman et al. (1995) zero-point for 
F814W propagates directly to systematic errors in the V; mag- 
nitudes. The F814W zero-point uncertainty is listed by Holtz- 
man et al. (1995) as ±0.006 mag. Their calibration, however, 
is explicitly valid for the spectral energy distributions of stars, 
so the calibration has larger uncertainties when applied to the 
SEDs of galaxies, probably on the order of 1-2%. Furthermore, 
the Holtzman et al. (1995) zero-point is valid for 0''5 radius 
aperture magnitudes. Our surface brightnesses will have been 
measured using data deconvolved with a point-spread function 
with a finite radius (1''5). Therefore, we must correct our data 
from an effective aperture of r = 1''5 to an effective aperture 
of r = 0''5. A reasonable estimate for the uncertainty in this 
correction is ±0.01 mag. Another source of error arises from 
uncertainties in the Extragalactic Dista nce Scale Key Project's 
correction to the zero-point (see, e.g., Stetson 1998); we esti- 
mate that it introduces an uncertainty on the order of ±0.025 
mag (half the effect). 

The color term in the transformation to redshifted V~ is small, 
with the result that errors in K due to uncertainties in the 
(F606W-F814W) colors are diminished by nearly a factor of 
five. These color errors include the uncertainty in the F606W 
zero-point, which we assume to be on the order of 0.05 mag, 
compared to the estimate by Holtzman et al. (1995) of ±0.005 
mag (the F606W zero-point was determined theoretically, and 
no additional calibration since Holtzman et al. (1995) has been 
attempted). With regard to the colors, any systematic varia- 
tions in the CCD gain and charge transfer, for example, cancel. 
As stated earlier, the error in the adopted Galactic reddening is 
likely to be small, no greater than ±0.01 mag. Thus, the pri- 
mary component of systematic error is the F606W zero-point. 
If the (F606W-F814W) color error can be thought of as the 
quadrature sum of the Holtzman et al. (1995) zero-point errors, 
then the resulting error in V- from galaxy colors is still :;0.01 
mag. Color gradients within the galaxies themselves could be 
a larger factor, but van Dokkum et al. (1998a) showed that, in 
the mean, the colors change very little as one adopts different 
standard radii for measurement of the colors. 

Errors in the coefficients in the transformation come from 



several sources. Uncertainties in the filter constants are esti- 
mated to be ±0.03 mag in the final transformation. Variations 
in the intrinsic galaxy SEDs lead to (slightly) different values 
for the color and constant terms in the transformations them- 
selves, introducing an uncertainty of ±0.01 mag. Lastly, adop- 
tion of the Hayes (1985) spectrum of Vega may also lead to 
errors if the WFPC2 photometric calibration sources were not 
strictly on the same spectrophotometric system as the Hayes 
spectrum. Unfortunately, the magnitude of this uncertainty is 
not well known, though we estimate it to be no larger than 
±0.05 mag based on the agreement of several sources of fil- 
ter constants with with Holtzman et al. (1995). 

By adding these contributions in quadrature, we estimate a 
total error of ±0.06 mag in V, and ±0.05 mag in (B--V,) color. 
For the analysis of the fundamental plane, these uncertainties 
can be treated as systematic errors. 

4. FITTING DE VAUCOULEURS PROFILES 

The photometric parameters, r^ and were obtained by 
fitting r'/'*-law profiles to the galaxies and two distinct ap- 
proaches were used to fit r'/'^-laws to the data. The first ap- 
proach involved restorat ion of the data (deconvolution) using 
the CLEAN procedure (hogbom 1974 ), isophote fitting using 
GALPHOT ( |j0rgensen et al. 1992 ), the derivation of a circu- 
larized, integrated surface brightness profile {i.e, a "growth" 
curve), and its fit with an integrated r'/'^-law. 

The second approach involved fitting two-dimensional 
gala xy models to the images directly (cf. v an Dokkum & Franx 
1996). The algorithm involves convolving two-dimensional 
r'/"*-laws with a PSF, iteratively adjusting the parameters until 
the model provides a good fit to the WFPC2 images. Nearby, 
contaminating objects can be handled elegantly, since multi- 
ple objects can be fit simultaneously (see also Simard 1998). 
For both of these procedures we used artificial PSFs for each 
g al a xy. These PSFs were created using Tiny Tim 4.0b ( Krist 
1995). We estimate the effects of PSF errors in §^ 

4. 1 . Fitting Integrated Surface Brightness Profiles 

In this section, we discuss the procedures used for fitting the 
integrated profiles derived from the WFPC2 images. Before de- 
riving the surface photometry, the images must be deconvolved 
because the structural parameters and can be seriously af- 
fected by the WFPC2 PSF and its broad wings. We chose to 
use the CLEAN met hod of deconvo lution because it is fast, ef- 
ficient, and simple ( |Hogbom 1974| ). The output images were 
convolved with a Gaussian smoothing function to reduce the 
noise in the output image. The residual maps typically con- 
tained ^ 1 % of the counts in the original image but were never- 
theless added to the output frames in order to preserve flux. We 
fitted elliptical is ophotes using the program GALPHOT (J 0r- 
gensen et al. 1992). The radii were spaced logarithmically, with 
an inner isophote radius of r = 0.3 pixels, extending out to radii 
typically five times larger than r^. 

In deriving structural parameters for our galaxies, we dupli- 
cated the procedures used on nearby galaxies, fitting the inte- 
grated profiles rather than the profiles themselves. 

4.1.1. Fitting for r^ 

The flux in the CLEANing process has been restored with an 
effective PSF, chosen to be a Gaussian with FWHM = 1 pixel. 
At any given radius, one needs to correct the analytical r'/"*- 
law growth curve for the convolution by the restoration PSF. 
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The enclosed flux within a given radius for an r^^'^-law profile 
is analytically simple: 



„7.676 



F(r) = nrih 



7(8,^) 



X = 1.616 



7.676-8 

27(8,^) 



7! 



(9) 



(10) 
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The least-squares fit can be linearized straightforwardly, since 
the derivatives of Equation |ll], with respect to and 1^ (or {I)e), 
are also analytically simple. 

For a given set of structural parameters (feje), one can read- 
ily compute how much interior flux should be redistributed out- 
side of a given radius in the growth curve. In this fashion, one 
adjusts the gr owth curve due to the smearing by the restoration 
PSF (see, e.g., [SagliaefaZ. 1993 ). At each radius, one computes 
a set of aperture corrections to correct the model to a "seeing"- 
convolved model (assuming circular symmetry). These correc- 
tions are determined numerically for each set of structural pa- 
rameters tested in the search for minimum. The expression 
for the convolved model growth curves is 



G(r) = 



2tt pI{p)P{k)Jo{kr) Ji {kp)dkdp (11) 



where the restoration PSF is described by its transform, P{k), 
and Jq and Ji are zeroth and first order Bessel functions. For 
I{p), we use the fitted profile. We use the Levenberg-Marquardt 
method to perform the minimization, fitting G(R) to the in- 
tegrated surface brightness profiles. Once the set of structural 
parameters has been found which minimizes x^, we have r^,, 
le, and {I)e for the fitted model. Once the model has been fit, 
the galaxy profile may not be perfectly fit by the model within 
re, the radius within which (7)^ is defined. Therefore, we use 
the observed i.e., the enclosed flux within r^, for the ob- 
served profile instead of the enclosed flux for the model, and 
have corrected this surface brightness for the smearing effects 
of the restoration PSF (see Equation O). Using either the ob- 
served or model (7)^ has an effect on r^I)^'^^ of less than 0.1% 
in the mean. 

In Figure |^, we show the calibrated F814W surface bright- 
ness profiles and integrated surface brightness profiles. The fit 
by the integrated de Vaucouleurs profiles are shown by thick 
solid lines in the right-hand panels. The implied r'/^^-law pro- 
files are shown in the left-hand panels. The residuals are shown 
in the lower panels, using thick solid lines as well. In general, 
the quality of the fit to the integrated profiles is good. 

4.1.2. 



PSF may not be perfect. These effects are exacerbated by the 
under-sampling of the data. 

The outer radius was chosen in order to minimize the effects 
of sky-subtraction errors. The shape of a galaxy's outer profile 
can be greatly affected by errors in the sky background, even 
though fitting of a parameterized growth curve is intended to be 
robust against such uncertainties. 

We find that the results are quite stable to the choice of in- 
ner and outer fitting radii. We varied the outer fitting radius 
between 2.0, 2.5 and 3rg, and the inner fitting radius between 
0, 0.5, 1 pixel, and 2 pixels. (For Hq = 65 kms-'Mpc"' and 
^0 = 0.1, one WFPC2 pixel corresponds to 0.48 kpc, equiva- 
lent to 1" at the distance of Coma.) For all combinations of the 
fitting range, mean and median differences in are less than 
±1.5%. We found that Irg was a suitable compromise for the 
outer radius in fitting the growth curve. After investigating the 
uncertainties due to PSF errors and the under-sampling (§|1), we 
adopted an inner fitting radius of r = 1 pixel. 

In the extreme case, changing the inner and outer fitting radii 
to r = and r = Sr^., the mean systematic offset in is about 
1%, with an rms scatter in the values of 9%. The system- 
atic offset in rj^ ''^, the combination of the structural parame- 
ters which appears in the fundamental plane, is negligible, with 
an rms scatter of 0.1%. The stability of the effective radii to 
the inner fitting boundary gives us confidence that the proce- 
dure outlined above properly corrects for the convolution by 
the restoration PSF. 

4. 1 .3. Sensitivity to Sky Subtraction 

One important source of error not addressed by varying the 
inner and outer fitting radii was uncertainties in the determi- 
nation of the sky background. When the imaging data were 
processed in van Dokkum et al. (1998a), the background was 
modeled and subtracted using order n ^ 15 polynomials. How- 
ever, for any given galaxy there may still be residual back- 
ground fluctuations that were not accurately removed by the 
modeling. Thus, we followed the procedure in J0rgensen et 
al. (1992), where low-order power laws (m = 2, 3) are fit to the 
outer parts of the profiles using an unweighted least-squares 
routine. The sky value is estimated from the asymptotic con- 
vergence of 7 = 7or-"' + sky. 

The rms scatter in extrapolated sky values was 5e- (-24 
mag/n"). By ignoring this fit for the local background, the ef- 
fect on r^ is about 1%, in the mean, with a standard deviation of 
3%. No effect on the fundamental plane parameters was found. 

In summary, the fundamental plane parameters derived from 
the fitting of the integrated profiles, for all combinations of the 
inner and outer fitting radii, are extremely robust against the 
choice of fitting range. The insensitivity to the fitting range 
and to uncertainties in the sky subtraction suggest to us that the 
growth-curve parameters are well-determined. 



Sensitivity to Fitting Range 

In the fitting process, one must select a range of radii over 
which to fit the integrated surface brightness profiles. At small 
radii, the profile shapes are affected by PSF errors which may 
have been exacerbated by the (de)convolution procedures. Sec- 
ondly, at large radii the profiles themselves become noisy, and 
thus the growth curve shapes can become suspect. 

Three effects limit the choice of inner radius in the fit: (1) 
the CLEAN procedure may not be perfect and noiseless; (2) 
the Tiny Tim construction of the WFPC2 PSF may not be per- 
fect; and (3) the correction for the smearing by the restoration 

*In order to linearize tlie fitting procedure, we continue to describe the growth curve's derivatives by their analytical expressions 
G(R)/F(R), where G and F are the convolved and unconvolved values of the cumulative flux. 



4.2. Fitting the Images Directly 

In this section, we derive structural parameters by fitting two- 
dimensional model galaxies d irectly to the WFPC2 images (see 
van Dokkum & Franx 199^ for a discussion of the technique 
used here, and oth e r authors for variants of the method e.g.. 



Schade et al. 1997 , Simard 1998 ). Images of model galaxies, 

multiplied by the ratio of 
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for a given set of structural parameters, are convolved with the 
Tiny Tim PSFs and compared directly to the data frames. The 
structural parameters are adjusted until x^, or some other mea- 
sure of the goodness of fit, is minimized. The local topology of 

is generally used to estimate the uncertainties in the fit. 

Procedures for fitting the images directly have the advantage 
that one performs a convolution rather than deconvolution. The 
advantages of such a process include the preservation of the 
image's noise characteristics, and the fitting of multiple objects 
simultaneously. The current disadvantage is that the 2D mod- 
els are built using a single, global ellipticity and position angle 
for each galaxy. Nonetheless, as the apparent information con- 
tent continues to shrink with the decreasing angular sizes of 
high redshift galaxies, useful structural parameters are still ex- 
tractable via such modeling, so long as the PSF is reasonably 
well known. 

Such fitting methods assume a functional form to the galaxy 
profile, such as an r''''*-law, but can be expanded to include mul- 
tiple components (e.g., r'/"*-law bulge plus exponential disk). 
The parameters of the least-squares fit to an r'/"*-law are: (1) 
(x,y) position of the center; (2) e, the apparent ellipticity; (3) 
position angle; (4) effective radius; (5) surface brightness at 
the effective radius; and (6) the sky background. Multiple- 
component (e.g., bulge-plus-disk) fits require additional param- 
eters and complicated minimization algorithms, but even such 
expanded profile shapes do not yet take into account devia- 
tions from purely elliptical isophote shapes, or twisting of the 
isophotes with radius. 

Using the fitting procedure of van Dokkum & Franx (1996), 
we derived the r'/'^-law structural parameters for the galaxies. 
Given the potential uncertainties in the processing of the sky 
background in the CCD images, we performed the fitting first 
with the sky background fixed at zero, and a second time while 
fitting for the local sky background. When the sky background 
is fit, the re changes by a mean (median) of +2% (+0.1%), with 
a scatter of nearly ±10%. The additional fitting of local sky has 
a negligible effect on the fundamental plane parameter in the 
mean (< 0.1%), with a scatter of 0.5%. For some galaxies, 
could change by as much as 30-40%, but even for these extreme 
cases the fundamental plane parameter only changes by ^ 2%. 

In Figure we show the galaxy images, after subtraction 
with the best-fit 2D r'/'^-law models. Many disk-like structures 
are clearly visible in these residuals and bulge-plus-disk fits to 
the ID profiles are discussed in §^ 

4.3. Comparison of Growth Curve and Image Modeling 

Ideally the two fitting procedures should yield similar, or 
even identical structural parameters. In practice this will not be 
the case for several reasons. For example, the two procedures 
have intrinsically different weighting schemes in the fitting pro- 
cess. Furthermore, galaxies do not have surface brightness pro- 
files which obey perfectly an r'/'^-law, or any other simply pa- 
rameterized function. While the internal tests listed above, such 
as varying the fitting range of the integrated surface brightness 
profiles, or allowing the image modeling to fit the sky back- 
ground, help estimate some of the uncertainties in the measure- 
ments, the methods themselves suffer from their own systemat- 
ics. Therefore, we now compare the structural parameters de- 
rived from the two methods in order to help estimate remaining 
systematic uncertainties. 

The comparison of effective radii derived by the two algo- 
rithms is surprisingly good, considering the range of apparent 
morphologies covered by the sample. The median difference 



in effective radii between the two algorithms is only 1%. The 
standard deviation is 13% in r^, precisely the quadrature sum of 
the rms scatter derived in the individual internal tests. The me- 
dian offset in effective surface brightness is -4%. Thus, there is 
a systematic offset in the effective fundamental plane parame- 
ter of 1 .5%, with a standard deviation of 2.5%. Nominally, one 
would expect that the error in the fundamental plane parame- 
ter should be smaller than the error in the individual structural 
parameters, because of the correlation between r^ and surface 
brightness. However, the two fitting methods weight the data 
differently and are therefore affected differently by the under- 
sampling and PSF errors, which are addressed in the next sec- 
tion. Nevertheless, these results indicate that the errors in the 
fundamental plane parameters are small. 

5. REMAINING INSTRUMENTAL EFFECTS 

While the above tests show that the fundamental plane pa- 
rameters appear well constrained, there are potential sources 
of error in common to both methods. The two most important 
effects are (1) the under-sampling; and (2) the use of model 
point-spread functions. Even those these effects are common to 
both methods, the data are weighted differently in the cal- 
culation. As a result, the two fitting methods are not affected 
identically. 

5.1. The Ejfects of Under- sampling 

In order to test the effects of under-sampling, and the depen- 
dence of the PSF on sub-pixel positioning, we generated a lOx 
subsampled Tiny Tim PSF. We created a set of noiseless sub- 
sampled r'/'^-law galaxy images using the subsampled PSFs. 
We then created two sets of WF CCD sampled versions of these 
model galaxies. The first set were simply block-averaged to 
the nominal WF CCD sampling and c onvolved with a charge- 
smearing kernel (Krist & Hook 1997). Thus, the galaxy cen- 
troids were located at the precise corner of a WF CCD pixel. 
The second set were created by shifting the subsampled model 
galaxy images by 0"05 (half a WF CCD pixel), then block- 
averaged, and convolved with the charge-smearing kernel. In 
each case, the block-averaged subsampled PSFs are equivalent 
to the (under-)sampled, real ones. Thus, we had two sets of 
model galaxy images which differed in position by half a WF 
CCD pixel. We then used the nominally sampled Tiny Tim 
PSF, centered at the integer pixel location, to derive structural 
parameters for these galaxies. 

The model images were analyzed in the same manner as the 
CL1358+62 galaxies, as described in §4.1 and §4.2. The sys- 
tematic errors in r^ varied with galaxy size, as did the errors in 
rJe ''^. For = 4 pixels, the error in the fundamental plane pa- 
rameter was only -2%. For r^, = 2 pixels, the error is -5%. For 
these tests, we restricted the growth curve fitting to radii r > 1 
pixel. By fitting to the very center, the systematic error in rj^-^^ 
was nearly -10% for galaxies as small as r^, = 2 pixels. Even by 
directly fitting the images with the two-dimensional models, the 
derived structural parameters are susceptible to errors due to the 
under-sampling. For galaxies with r^ = 1 pixel, the systematic 
errors due to under-sampling can be as large as -50% in r^, and 
^ 10% in rglg '^^. Fortunately, the smallest galaxy in our sample 
has r^, w 2 pixels; the median effective radius is « 6 pixels. 
Therefore, we anticipate that the net systematic effect of the 
under-sampling on the fundamental plane is small, being 1-2% 
on average. 

We caution against fitting profiles or growth curves too close 
to the centers of galaxies, even when deconvolution of the 
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WFPC2 images has taken place, because the under-sampling 
makes fitting small galaxies quite problematic. Techniques 
which impro ve the sampling, such as "drizzling" (H ook & 
Fruchter 1997) should improve the usefulness of the data for 
such objects. 

5.2. A Comparison of Stellar and Model Point-Spread 
Functions 

In both procedures, we have reli ed on the a ssumption that 
the Tiny Tim point-spread function ( Krist 1995 ) is an accurate 
match to the PSF of the observations. Even when this assump- 
tion is be correct for a single exposure, the process of combin- 
ing multiple exposures in each pointing introduces extra "jit- 
ter" in the PSF of the resulting image. The situation is actually 
worse because one has only a finite number of images, e.g., 
N = 3, and not the limit of — + cx3 images. Therefore the indi- 
vidual jitter vectors are not going to be isotropically distributed, 
for each galaxy (or star) in a given pointing. 

Therefore, we tested to what extent our results are effected 
by errors in the model PSF. We built model galaxy images us- 
ing a bright star as the observed PSF, and then analyzed these 
images using a Tiny Tim model PSF. The star was located at 
approximately (550,541) in CCD #4 of pointing #1. We per- 
formed this test using observed stellar and Tiny Tim PSFs with 
1"5 diameter 

We generated a lOx subsampled model PSF at the integral 
pixel location of the star, shifted and resampled it to the pixel 
size of the WF CCDs, and then convolved with the charge- 
smearing kernel. We measured the flux-weighted centroid of 
this nominally sampled PSF and adjusted the sub-pixel shift 
until its centroid coincided with that of the star. The effects due 
to under-sampling in the star and Tiny Tim point-spread func- 
tions should cancel out, because they are positioned at identical 
sub-pixel locations (to within 5% of a pixel). 

We created a set of 2D model galaxy images which had been 
convolved with the stellar image. The models had pure r'/^- 
law profiles with effective radii G (1,2,4,8, 16) pixels. Us- 
ing the sub-pixel shifted Tiny Tim PSF, we then measured the 
structural parameters of the model using the direct image fit- 
ting. The measured effective radii underestimated the correct 
values by 17%, 12%, 8%, 6%, and 4%, respectively. This sug- 
gests that the model PSF was broader than the star, but because 
of the under-sampling, this difference in width was not initially 
obvious. The fundamental plane parameters were systemati- 
cally larger by 6.8%, 3.6%, 1.8%, 0.8%, and 0.2%, respectively. 
These offsets are independent of model ellipticity. 

The smallest galaxy in our sample has an effective radius of 
r<; « 2 pixels (0"2). The median half-light radius of the sam- 
ple is 0''6, with 25 and 75 percentiles of 0".35 and ~1". Thus, 
for the upper half of the sample, the uncertainties on the fun- 
damental plane parameters, due to improper matching of the 
point-spread functions, are < 1%. For the 25% of the sample 
with 0''35 < rg < 0".6, the effects are 3%. For the remainder, 
the uncertainty is < 4%. For galaxies at higher redshift, or for 
galaxies farther down the luminosity function, one must take 
extreme care to match the point-spread function to the data, 
given the potentially large systematic uncertainties at r^, w 1 
pixels. 

Using two distinct fitting methods, we have reliably derived 
effective radii and surface brightnesses. The uncertainties in 
the effective radii of the de Vaucouleurs profile fits are small, 
at the level of -13% rms, with a mean offset of only 1%. The 



measurements are robust against fitting range, uncertainties in 
the sky background, and fitting method. The uncertainties in 
the fundamental plane parameters themselves is even smaller, 
-2.5% rms, with potential systematic uncertainties at the level 
of < 3%, due to errors in the PSF, under-sampling, and differ- 
ences between the fitting methods. 

However, other sources of uncertainty in half-light radii may 
persist, due to deviations in galaxy profiles from the r'/"*-law 
form. Therefore, we continue with a discussion of systematic 
deviations from the r'/'^-law. 

6. SERSIC R'/'^-LAW profiles 

The actual surface brightness profiles of galaxies often de- 
viate from pure r'/'^-laws. In the sample presented here, both 
bulge and disk components are readily visible in many of the 
images (see Figures |l| and We now investigate the fitting 
of more generalized surface brightness profiles to our data to 
test (1) if a more complex profile is better suited to the pro- 
files of the CL1358H-62 galaxies; (2) if such departures from 
an r'/'^-law are important for the determination of their struc- 
tural parameters; and (3) if they allow for accurate and objective 
morphological classification. 

In this section, we generalize the parameterization of the 
surface brightn ess profiles to the r'/"-law (C aon, Capaccioli, 
D'onofrio 1993, praham & Colless iwf^ where n, in the case 
of this paper, is restricted to the set of integers n G {1 . . .6}. 
Note that the case of n = 1 corresponds to an exponential disk, 
and n = 4 is equal to the de Vaucouleurs profile. 

The inclusion of n serves two purposes. First, the best-fit n is 
a a quantitative descriptor of galaxy morphology. Graham et al. 
(1997) found that the giant cluster ellipticals and BCGs follow 
n > 4 profiles. Furthermore, Saglia et al. (1997) showed that 
r'/" is equivalent to a subset of bulge-plus-disk models, though 
the relationship between n and B/D is not a straightforward 
one. 

6. 1 . Fitting the Growth Curves 

Integrated r'/"-laws, wit h n = 1,2,3,4.5,6 have been fit to 
the growth curves of §4.1. The enclosed flux within a given 
radius for an r'/"-law profile can be simply expressed: 



F„(r) = nrl{I)„ 



2"f(2n,x) 
(2«-l)! ' 



(12) 



x = i2n-Q.324)[ - 

. r„ 



i/« 



where r„ refers to the half-light radius of the corresponding 
r^/"-law fit. Fundamental plane parameters will be expressed 
as r„I^ '^^ andr„(/}0 ™. 

As in the case of the de Vaucouleurs fitting (n = 4), we adopt 
structural parameters derived from an inner fitting radius of 1 
pixel, and an outer radius of 2r„. We tested the structural pa- 
rameters from these fitting functions with direct image model 
fitting using the r'''"-laws for a restricted set of n (n = 1 . . .4). 
Comparisons of the structural parameters derived from the im- 
age and growth-curve fitting for n = 1,2 , 3 were similar to the 



comparisons for n = 4 discussed in §4.3 (in which differences 



in r4l^ ''^ were 1.5% in the mean, with 2.5% rms scatter). 

6.2. A Comparison of the r'/" Profiles 

By expanding our growth curve fitting to a range of pro- 
file shapes, one expects to gain additional knowledge about the 
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galaxy sizes and morphologies. In this section, we investigate 
whether or not the new measures of half-light radii are necessar- 
ily more accurate than those obtained with the de Vaucouleurs 
profile. Figure ^a) shows the distribution of the ratio r„/rj, 
(where = rn) using the best-fit n for a given galaxy. This 
distribution is quite broad, suggesting that the fit of a pure de 
Vaucouleurs profile can lead to large errors in half-light radius, 
similar to what has been found for nearby samples (C aon, Ca- 
paccioli, D'onofrio 1993). 

For a given n, the ratio of /V4 appears to be a well defined 
function. As shown in Figure a), ln(r„/r4) « -1.32 + 0.33«, 
almost irrespective of the true shape of the galaxies, with a scat- 
ter of less -5%. Thus, the fit of an n = 6 profile to an intrin- 
sically n = 2 galaxy will find too large a half-light radius by 
nearly a factor of three. As a result, a large bulk of an n = 6 
profile is fit to parts of the growth curve which are more sus- 
ceptible to errors in the sky subtraction. In crowded regions, 
where a large source of systematic error is the background de- 
termination, residual background light can easily appear as a 
shallow envelope in the integrated surface brightness profiles, 
masquerading as the profile of an r'/*-law, or worse yet, even 
larger values of n. Of course, the other extreme is also possible. 
By over-estimating the sky, one can effectively decrease the n 
derived for a given galaxy. Therefore, when expanding the pro- 
file fitting to include the overall shape («), one must treat the 
background subtraction with care. 

While for fixed values of n, the individual parameters of r„ 
and /„ are not particularly sensitive to errors in background sub- 
traction (see the earlier case of n = 4), the chief uncertainty 
is which value of n to assign to a given galaxy. By varying 
the range of radii over which one compares the different r'/" 
fits, one begins to test the sensitivity of n to the uncertainties in 
sky subtraction. A 3% error in sky subtraction (§^ can lead to 
changes in n. An sa 1, and errors in half-light radius of 20-30%. 
The typical scatter introduced in n due to the sky extrapolation 
alone is of order A« « ±0.2. Nearly a third of the sample, 
however, have An values larger than this, and many of these 
are the « * 4 galaxies for which uncertainties in n on the order 
of unity are not uncommon. We adopt An as the uncertainty 
in the measurement of n, and the faintest galaxies, as expected, 
tend to have the largest uncertainties. We conclude that the r'/"- 
law does not yield more accurate estimates of half-light radii, 
especially for galaxies fit by n > 4. 

Several authors have commented upon a trend of galaxy 
shape paramet er, "n," with galaxy size (C aon, Capaccioli, 
D'onofrio 1993, praham & CoUess 1997 ). However, in Fig- 
ure ^(a) we saw that the ratio of r„/r4 scales in a well-defined 
manner By fitting a low-n profile to a galaxy, one will automat- 
ically find a smaller half-light radius than when fitting a larger 
n profile. The half-light radii from low-n fits are simply smaller 
than the half-light radii one derives when fitting larger-n pro- 
files, independent of the true shape of the galaxy profile. This 
bias may complicate the nature of the trend of n with galaxy 
size. 

6.3. Effective and Mean Effective Surface Brightness 

Ideally, the fundamental plane should be insensitive to the 
detailed profiles of (early-type) galaxies. However, the funda- 
mental plane parameter constructed with the surface brightness 
at an effective radius is not robust. As can be seen from Fig- 
ure |(b), the ratio of ri/j'^*' to re/^-^^ is a factor of -2.5. The 
ratios of the r„lf^-''^ / r^I^^^^ to each other are independent of the 
real shapes of the galaxy profiles, to better than a few percent. 



In Figure ^c), however, we see that the fundamental plane pa- 
rameter constructed with the mean surface brightness (within 
r„) is much more robust with respect to n, showing only -10% 
peak-to-peak difference in r„(/)5]^*. The typical scatter about 
this relation is only a few percent. 

The ratio of r„ (Z),',''^^ to r^ (Z)"'^^ is virtually independent of 
galaxy profile as well. This indicates that every n = 2 galaxy fit 
by a de Vaucouleurs profile is actually measured incorrectly in 
precisely the same way as every other n = 2 galaxy, independent 
of galaxy size. 

In summary, the fitting of the Sersic r'/"-law profiles does not 
lead to more accurate photometric parameters for the purposes 
of fundamental plane analysis. The derived half-light radii are 
correlated with the best-fit n, which is sensitive to errors in the 
sky background. Fortunately, the fundamental plane is insen- 
sitive to these uncertainties in the photometric parameters, be- 
cause the fundamental plane parameters, when constructed us- 
ing mean effective surface brightnesses, remain constant at the 
level of 3%. 

7. BULGE/DISK DECOMPOSITIONS AND QUANTITATIVE 
MEASURES OF MORPHOLOGY 

Classifying galaxies by morphology becomes increasingly 
difficult at high redshifts. Distant galaxies are typed using 
images taken with (linear) CCD detectors, while nearby sam- 
ples have been primarily typed using (non-linear) photographic 
plates. Furthermore, the distant galaxies typically span fewer 
resolution elements, and the images have poorer S/N than im- 
ages of the nearby prototype elliptical and lenticular galaxies 
with which we wish to make comparisons. In order to maxi- 
mize the use of our data in determining galaxy morphologies, 
we classify the galaxies in our sample by both visual and quan- 
titative morphology. The visual classifications have been taken 
from Fabricant et al. (1999). 

In the previous section, we derived the n parameters for the 
galaxies in our sample. As stated earlier, the Sersic (1968) 
n profiles are a "r estricted" subset of bulge-plus-disk models 
(Saglia et al. 1997). In fitting the growth curves, the use of n as 
a descriptor of morphology is formally more robust than bulge- 
disk decompositions. One needs only three quantities to de- 
scribe a profile, (n, r^,/;,), rather than four, ('^/2'^f/2''f/2'^i/2)- 

In this section, we test whether n reasonably reflects B/D ra- 
tios for our sample of galaxies. To this end, two-component, 
bulge-plus-disk growth curves were fit to the integrated surface 
brightness profiles of our galaxies, using the same fitting radii 
as above. The growth curves were fit between 1 pixel (0"1) 
and 2r„, defined by the "best" fit r'/"-law (above). As in the 
fitting of the r'/'^-law growth curves, we used the Levenberg- 
Marquardt method to search for the parameters which mini- 
mized x^, taking extra care to avoid local minima. 

In Figure |^, we show the fitted bulge-plus-disk superposi- 
tions and the residuals from this fit using thin solid lines. As 
is shown in the figures, the galaxies in our sample appear to be 
well described by the superpositions of de Vaucouleurs bulge 
(n = 4) and exponential disk (n = 1) profiles. The individual 
bulge and disk components are plotted as dashed and dotted 
lines, respectively. Note that the bulge-plus-disk models tend 
to provide better fits than a pure r'/'^-law. The addition of an 
exponential disk improves the fit for many of the galaxies. The 
Eh-A galaxies, ID # 209, 328, and 343, are bulge-dominated, 
but clearly have disks at the level of 10-35%. Over the en- 
tire sample, the median bulge fraction is 84%. Upon visual in- 
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spection of Figure ^ one finds disk-like structures in about half 
of the sample, consistent with the visual classifications. Spiral 
arms are also plainly visible in many, including two of the E+A 
galaxies. 

In many of the galaxies, there appears to be evidence for cen- 
tral exponential (disk) structures. Most of the galaxies with 
larger bulge fractions show evidence for central exponential 
components in their surface brightness profiles. This does not 
appear to be an artifact of the fitting procedures or errors in 
the PSF, as many of these central "disk" components persist at 
radii r^l pixel. Evidence of such structures might appear in 
the kinematic profiles but seeing would c omplicate the inter- 
pretation of such data (Kelson et al. 1999a). A full study of the 
entire members hip catalog using two-dimensional model fitting 
of the images (Tran et al. 2000), would help confirm whether 
this abundance of central exponentials is the norm, or an artifact 
of the deconvolution process. 

Table |l] lists the optical morphologies (Fabricant et al. 1999) 
with the best-fit n parameters. We list (1) galaxy ID, (2) vi- 
sual morphology, T, from Fabricant et al. (1999), (3) (B-V) 
color, also from van Dokkum et al. (1998a), (4) r^, (5) {n)e, 
(6) r,(/}0■^^ (7) «, (8) r«/2, (9) {fiYl^^, and (10) M'''. The 
fundamental plane parameters have been expressed using units 
of kpc and Lq/pc^ (Hq = 65 kms"'Mpc"', = 0.1). 

Table | gives the best-fit bulge fractions for the galaxies in 
this sample, and lists the individual values of (3) ri/2, (4) {^j)\i2, 
(5) ri/2{t)lj2- (6) ri/2 for the bulge component, (7) (/i}i/2 for 
the bulge component, (8) /2 for the disk component, (9) (fj.) i /2 
for the disk component. Bulge Fraction within (10) r < Q".2, 
(11) r < 0"5, and (12) r < Ir^/i, (13) total bulge fraction, and 
(14) luminosity-weighted ellipticity, are given for each galaxy. 
Scale lengths are listed in arcsec, and surface brightnesses in 
V; mag/n". The tables include only formal errors and do not 
include other sources of errors, such as uncertainties in the cal- 
ibration or errors due to PSF mismatch. 

There were two explicit reasons for expanding the form of 
the fitted surface brightness profiles to the r'''"-laws and bulge- 
plus-disk superpositions; they are to (1) determine accurate 
measures of half-light radii; and (2) objective morphological 
classifications. Earlier we compared the half-light radii from 
the Sersic profiles to those derived using the de Vaucouleurs 
profile. In Figure ^b), we show the distribution of bulge-plus- 
disk half-light radius, rg+o, divided by r,., the half-light radius 
of the de Vaucouleurs profile fit. The distribution is broad, and 
with Figure H(a), should further caution the reader about the 
reliability of half-light radii in general. 

The other reason behind expanding the model surface bright- 
ness profiles to Sersic profiles, and bulge-plus-disk superposi- 
tions was to test whether these forms could be used to objec- 
tively classify the galaxies, given the recent interest in quanti- 
tative morphologies. In Figure ^ we compare visual and quan- 
titative estimates of morphology to help test their reliability. 
Figures 0(a-c) compare BF, the bulge fraction of total light, 
n, the best-fit r'''"-law, and e, the flux-weighted apparent ellip- 
ticities from the image modeling to the visual classifications, 
T, of Fabricant et al. (1999). In general, the visually classified 
E/SOs (T < -3) have high bulge fractions. Those galaxies with 
n > 4 are all of type SO/a and earlier. While T and n are clearly 
correlated, the relationship between profile shape, n, and visual 
morphology has large scatter. Figure ^d,e) shows BF and n 
as functions of e. In Figure ^f) one sees that galaxies with 
n * 4 typically have high bulge fractions. Clearly the relation- 



ship between bulge fraction, n, and optical morphology is not 
straightforward. That there appears to be no tight relation be- 
tween n and BF, however, should not be surprising. Saglia et al. 
(1997) showed that n is a complicated combination of the bulge 
fraction and the ratio of disk to bulge half-light radii. We con- 
clude that there is only a level of correspondence between the 
different estimates of morphology which allows one to discrim- 
inate between bulge-dominated and disk-dominated galaxies. It 
is clear from the figures that these quantitative measures are not 
successful at distinguishing subtle differences between adjacent 
morphological classes (e.g., Es from SOs). 

8. SUMMARY OF THE STRUCTURAL PARAMETERS 

The primary goal of this paper is to provide the photomet- 
ric parameters required to compare the fundamental plane in 
CL1358H-62 to that measured locally. We have shown that the 
growth-curve structural parameters are consistent with those 
derived from directly fitting 2D models to the images. 

However, the direct image fitting implemented by van 
Dokkum & Franx (1996) is currently restricted to rather simple 
galaxy models: a single component with constant ellipticity and 
position angle. Other, more sophisticated algorithms, such as 
those described by Simard et al. (1997), could be used as well. 
Given the magnitudes of the uncertainties in half-light radii (cf. 
Figure ||), the more sophisticated methods will not likely result 
in more accurate half-light radii, but simply in greater accuracy 
in the estimation of their errors. More distant galaxies, with 
smaller angular sizes, will require more sophisticated modeling 
to avoid serious systematic errors in structural parameters due 
to the under-sampling. Despite these issues, we have shown 
that the measurements of are quite robust. 

The current errors in the structural parameters, as deduced by 
the differences in the results of our two methods, give us confi- 
dence that more complicated procedures would not greatly im- 
prove on the analysis. What we do find, however, is that one can 
use n to grossly distinguish between bulge- and disk-dominated 
systems. These r'''"-laws, however, do not give much insight 
into any particular physical properties of the galaxies, such as 
ratios of bulge to disk radii, or bulge fraction. 

Given the uncertainties in the actual profile shapes of the 
CL1358H-62 galaxies, the fundamental plane parameter must be 
constructed with care. Use of the surface brightness at the half- 
light radius appears to be a strong function of the profile chosen 
for the fitting (see Figure |^). Adopting the mean surface bright- 
ness within the half-light radius produces a fundamental plane 
parameter which is robust against the adopted profile shape. 

Although we have outlined potentially large uncertainties in 
half-light radii, the correlation between r and (/) ensures that 
the fundamental plane parameter, r{I)^-''^ is conserved, with a 
scatter of 3%. In Figure ^ we show the correlation between 
(/) and r explicitly for the entire sample. For every trial fitting 
range and r'''"-law, we plot the mean surface brightness within 
the fitted half-light radius, normalized by the adopted mean ef- 
fective surface brightness from the de Vaucouleurs profile fit, 
versus galaxy radius, normalized by the effective radius from 
the de Vaucouleurs profile fit. The data do not deviate signifi- 
cantly from the correlation expected from pure r'/"*-law growth 
curves, shown by the dashed line. The slope running parallel 
to the fundamental plane, r{I)^-''^ = constant, is shown by the 
solid line. Given the large uncertainties in half-light radii, the 
strong error correlation may not only be coincidentally paral- 
lel to the fundamental plane, but also hkely seriously biases the 
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meas ured slope of the scaling relation (iKelson 1998L Kelson 
1999). 

Because the analysis most closely parallels earlier, ground- 
based work, we will use the results from the growth curve fit- 
ting in th e derivation of the fundamental plane (see K elson et 
al. 1999b). 

9. CONCLUSIONS 

Using deep HST imaging, we have extracted growth curves 
for a sample of 53 galaxies in the cluster CL1358H-62 at z = 
0.33. These growth curves have been fit using integrated r^/"'- 
laws, integrated Sersic r'/"-law profiles, and bulge-plus-disk 
superpositions. Comparisons of the fits between the different 
parameterizations of the profiles show an nns scatter of 30-40% 
in half-light radii, although the formal uncertainties in the half- 
light radii from a given fit are typically smaller by factors of 
three. The half-light radii derived from the Sersic profiles agree 
in the median with those found using the r'/'^-laws at the level 
of 1 %. The half-light radii derived using the bulge-plus-disk su- 
perpositions are systematically larger by 10%. Despite the po- 
tentially large uncertainties in the half-light radii, the quantity 
that enters the fundamental plane, r„{I)^^''^ has been shown to 
be quite stable. This combination is shown to be nearly parallel 
to the correlation between radius and mean surface brightness 
for an r'/'^-law growth curve. 

The median half-light radius of this sample is 0''61 (2.9 kpc, 
Hq = 65kms"'Mpc"', = 0.1), with a range of 1.26 dex in r^. 
(Using B/D decompositions, or r'/"-laws, only changes the me- 
dian by ±0''01.) The range is smaller than the 1.6 dex spanned 
by the J0rgensen et al. (1993) Coma sample, probably due to 
differences in selection criteria. 

The median bulge fraction is about 84%. Many of the galax- 
ies with larger bulge fractions show evidence for central ex- 
ponential components in their surface brightness profiles. The 



few E-hA galaxies in this sample invariably have disks (see also 
Franx et al. 1998. Their bulge f ractions are -65-90% (see also, 
e.g., Wirth, Koo, & Kron 1994). Half of the sample shows ev- 
idence for disky structures after subtracting two-dimensional 
r'/'^-law models from the galaxy images. This result is consis- 
tent with the visual classifications in which nearly half of the 
sample are of type SO/a trough Sb. Two of the three Eh-A galax- 
ies have evidence of spiral arms in their residual maps. 

In the future, we plan to derive accurate structural param- 
eters for the entire sample of confirmed cluster members and 
investigate the global correlations of these parameters with the 
signatures of stellar populations (col or), and environm ent (clus- 
tercentric radius and local densities) (Tran et al. 2000). The im- 



portance of these issues for the fundamental plane is discussed 
by Kelson et al. (1999b). Studies of galactic structural param- 
eters and scaling relations in distant clusters, like this program, 
clearly require wide-field imaging of the distant cluster fields. 
Two more WFPC2 mosaics have now been taken, of MS2053- 
04 (z = 0.58) and MS 1054-03 (z = 0.83). With these clusters, 
the evolution of the properties of cluster galaxies can be accu- 
rately mapped to redshifts approaching unity and the Advanced 
Camera for Surveys is expected to extend such studies to red- 
shifts beyond z = 1 . 
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Table 1 

ID r^l^- AND r^/"-LAw Structural Parameters 



















1 /n 






ID 


Type 


(S-V)z 


re 














(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


095 


SO 


0.894 


0.356 
±0.004 


19.84 
±0.02 


2.59 


5.2 
±1.2 


0.565 
±0.005 


20.59 
±0.02 


2.56 


110 


so/a 


0.850 


0.648 
±0.012 


20.32 
±0.04 


2.70 


2.3 
±0.4 


0.375 
±0.005 


19.53 
±0.03 


2.71 


129 


so/a 


0.830 


0.324 
±0.004 


19.27 
±0.02 


2.72 


4.2 
±0.3 


0.326 
±0.004 


19.28 
±0.03 


2.72 


135 


SO 


0.827 


0.351 
±0.009 


20.19 
±0.05 


2.48 


5.2 
±0.3 


0.493 
±0.010 


20.71 
±0.04 


2.47 


142 


so/a 


0.856 


1.414 

±0.046 


22.15 
±0.07 


2.49 


3.9 
±0.9 


1.398 
±0.046 


22.13 
±0.07 


2.49 


164 


so/a 


0.870 


1.165 
±0.035 


21.56 
±0.06 


2.58 


5.7 
±0.7 


1.955 
±0.100 


22.37 
±0.08 


2.56 


182 


SO 


0.839 


0.418 
±0.008 


20.70 
±0.04 


2.40 


3.2 
±0.3 


0.310 
±0.005 


20.25 
±0.04 


2.40 


209 


Sa 


0.601 


0.980 
±0018 


21.34 
±0.04 


2.57 


2.3 
±0.3 


0.557 
±0.005 


20.49 
±0.02 


259 


211 


SO 


0.870 


0.612 
±0.015 


20.58 
±0.05 


2.60 


3.3 
±0.4 


0.470 
±0.009 


20.18 
±0.04 


2.61 


212 


E 


0.868 


0.547 
±0.009 


20.44 

±0.03 


2.59 


5.2 
±0.3 


0.781 
±0.008 


21.01 

±0.02 


2.58 


215 


SO 


0.843 


0.338 
±0.006 


20.17 
±0.03 


2.47 


5.0 
±1.0 


0.460 
±0,007 


20.67 
±0.03 


2.45 


233 


E/SO 


0.847 


0.590 
±0.012 


19.95 
±0.04 


2.78 


5.3 
±0.4 


0.873 
±0,010 


20.58 
±0.03 


2.75 


234 


Sb 


0.689 


3.636 
±0.349 


23.68 
±0.18 


2.43 


3.0 
±1.0 


2.264 
±0.166 


22.97 
±0.16 


2.44 


236 


SO 


0.895 


0.539 
±0.009 


20.51 
±0.03 


2.57 


5.5 
±1.5 


0.892 
±0.010 


21.33 
±0.03 


2.53 


242 


E 


0.882 


1.201 

±0.026 


21.54 

±0.04 


2.60 


4.2 

±0.3 


1.210 
±0.026 


21.55 
±0.04 


2.60 


256 


E 


0.875 


0.956 
±0.010 


20.30 
±0.02 


2.88 


5.4 
±0.4 


1.542 
±0.010 


21.06 
±0.02 


2.86 


269 


E/SO 


0.913 


0.776 

±0,008 


20.05 
±0.02 


2.86 


5.0 
±1.0 


1 .056 
±0,008 


20.55 
±0.01 


2.85 


292 


so/a 


0.776 


0.911 
±0.031 


21.21 
±0.07 


2.58 


2.4 
±0.5 


0.562 
±0.009 


20.50 
±0.04 


2.59 


298 


SO 


0.914 


0.642 
±0.007 


19.93 
±0.02 


2.82 


2.5 
±0.5 


0.406 
±0.002 


19,23 
±0.01 


2.83 


300 


SO 


0.879 


0.386 
±0.004 


19.59 
±0.02 


2.70 


3.4 
±0.4 


0.310 
±0.002 


19.25 
±0.02 


2.71 


303 


E 


0.858 


0.494 
±0.007 


20.18 
±0.03 


2.63 


5.3 
±0.4 


0.729 
±0.007 


20.81 
±0.02 


2.61 


309 


WSO 


0.900 


0.433 

±0,003 


19.80 
±0.01 


2.69 


4.9 
±1.0 


0.572 

±0,003 


20.23 
±0.01 


2.68 


328 


Sa 


0.629 


0.878 
±0,016 


21.15 
±0.04 


2.58 


2.5 
±0.6 


0.550 
±0,008 


20.45 
±0.02 


2.59 


335 


so/a 


0.859 


0.512 
±0.009 


20.79 
±0.04 


2.46 


4.4 
±0.4 


0.551 
±0,009 


20.89 
±0.04 


2.46 


343 


SO 


0.670 


0.438 
±0.010 


20.66 
±0.05 


2.43 


3.3 
±0.4 


0.353 
±0.007 


20.33 
±0.05 


2.44 


353 


E^SO 


0.889 


0.958 
±0.007 


20.79 
±0.01 


2.73 


5.4 
±0.5 


1.495 
±0.008 


21.52 
±0.01 


2.70 


356 


Sa 


0.889 


0.815 
±0.011 


20.62 
±0.03 


2.71 


3.9 
±0.2 


0.741 
±0.011 


20.48 
±0.02 


2.71 


359 


SO 


0.851 


0.475 
±0.005 


19.95 
±0.02 


2.68 


4.4 

±0.4 


0.527 
±0,005 


20.12 

±0.02 


2.68 


360 


E 


0.861 


0.241 
±0,004 


19.65 
±0.03 


2.48 


5.5 
±0.6 


0.373 
±0,009 


20.37 
±0.04 


2.45 


366 


so/a 


0.863 


0.587 
±0,008 


20.68 
±0.03 


2.55 


3.3 
±0.4 


0.448 
±0,005 


20.28 
±0.03 


2.56 


368 


Sab 


0.862 


1.055 
±0.032 


22.09 
±0.06 


2.38 


2.4 
±2.6 


0.649 
±0.010 


21.34 
±0.04 


2.39 


369 


so/a 


0.879 


0.935 
±0.023 


21.70 
±0.05 


2.44 


4.3 
±0.4 


0.985 
±0.023 


21.78 
±0.05 


2.44 


371 


Sa 


0.871 


1.699 
±0.033 


21.84 

±0.04 


2.66 


3.5 
±1.5 


1.456 
±0.018 


21.61 
±0.03 


2.66 


372 


Sa 


0.867 


1.595 
±0,064 


22.09 
±0.08 


2.56 


3.2 
±1.2 


1.191 

±0,034 


21.64 
±0.06 


2.57 


375 


E 


0.894 


3.813 
±0,075 


22.44 
±0.04 


2.83 


6.0 
±0.2 


7.498 
±0,231 


23.53 
±0.06 


2.79 


381 


wso 


0.866 


0.392 
±0.008 


19.92 
±0.04 


2.61 


6.0 
±0.2 


0.756 
±0.020 


20.97 
±0.05 


2.57 


397 


so/a 


0.851 


0.727 
±0.020 


21.29 
±0.06 


2.46 


3.3 
±0.4 


0.472 
±0.011 


20.71 
±0.04 


2.45 


408 


SO 


0.861 


0.353 
±0.003 


19.70 
±0.02 


2.63 


3.4 
±0.5 


0.291 
±0.002 


19.41 
±0.02 


2.63 


409 


E 


0.895 


0.536 
±0.005 


21.38 
±0.02 


2.30 


4.1 
±1.1 


0.525 
±0.005 


21.35 
±0.02 


2.30 


410 


SO 


0.870 


0.497 

±0,008 


20.70 
±0.03 


2.47 


3.2 
±0.3 


0.370 

±0,005 


20.26 
±0.03 


2.48 


412 


E 


0.843 


0.368 
±0.007 


20.22 
±0.04 


2.49 


6.0 
±0.2 


0.795 
±0,023 


21.44 
±0.06 


2.45 


433 


so/a 


0.851 


0.371 
±0.011 


19.95 
±0.06 


2.57 


6.0 
±0.2 


0.907 
±0.039 


21.42 
±0.09 


2.51 


440 


so/a 


0.843 


0.663 


21.67 


2.30 


4.9 


0.843 


22.02 


2.30 


454 


so/a 


0.807 


0.983 
±0.014 


21.16 
±0.03 


2.63 


2.3 
±0.3 


0.613 
±0.005 


20.46 
±0.02 


2.64 


463 


SO 


0.883 


0.667 
±0.007 


20.50 
±0.02 


2.66 


3.2 
±0.3 


0.492 
±0.004 


20.03 
±0.02 


2.67 


465 


Sa 


0.873 


±0,009 


±0.03 




±0.3 


0.599 

±0,009 


20.98 
±0.03 


2.47 


481 


SO 


0.810 


0.478 
±0,015 


21.06 
±0.06 


2.35 


2.3 
±0.4 


0.334 
±0,005 


20.54 
±0.03 


2.35 


493 


E/SO 


0.809 


0.207 
±0,010 


20.09 
±0.10 


2.28 


4.1 
±1.1 


0.214 

±0,010 


20.15 
±0.10 


2.27 


523 


SO/a 


0.905 


0.910 
±0.016 


21.42 
±0.04 


2.52 


4.0 
±1.0 


0.863 
±0.016 


21.34 
±0.03 


2.52 


531 


E 


0.887 


1.048 
±0.009 


20.63 
±0.02 


2.82 


5.4 
±0.4 


1.682 
±0.010 


21.39 
±0.01 


2.79 


534 


E 


0.888 


0.464 

±0.008 


20.77 
±0.04 


2.42 


4.4 

±0.4 


0.486 
±0.008 


20.84 

±0.04 


2.42 


536 


E 


0.874 


0.890 
±0.009 


20.63 
±0.02 


2.75 


5.2 
±0.3 


1.355 
±0.010 


21.29 
±0.02 


2.73 


549 


Sab 


0.760 


1.592 
±0.058 


22.77 
±0.07 


2.35 


3.5 
±0.5 


1.380 
±0.033 


22.54 
±0.06 


2.36 


Note. 


— Morphologii-Lil 


V iLikcii Iroiii hibi 


cam ((/. I 1 


1L.'VU1-LI1I1L- iB-V 




lakcii iVoni vai 


D«kku]ii ( 


((/. (l'-)9,Sa)-, radii 


expressed 


in arcsec 


and surface 


brighlncsscs in m 


g per square 


arcsec, csccpl lor ihc 


undamc 


ilal plane parameters in which ihc radii have 



been transfonned to kpc and the surface brightnesses to Lq /pc ; these structural parameters were derived by fitting the int^rated profiles 
between r = o".l and r 2r^^2- 



Kelson et al. 



Table 2 

ID Bulge-Disk Decompositions 











T 1 


76 
'1/2 ^^>l/2 








D' k 




B 1 
u gi; 


F ■ l' 
VMi ion 






ID 


Type 


'-1/2 


(/^)l/2 


'■1/2 




'■1/2 


<m)i/2 


,, 

< ,2 


< .5 


< 2 '-1/2 


Toial 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


095 


SO 


0.704 


20.94 


2.55 


0.986 


21,84 


0,064 


17,75 


0,46 


0,64 


0.77 


0,85 


0,15 






±0.056 


±0.48 




±0.078 


±0.15 


±0,005 


±0,12 


±0,04 


±0,02 


±0,01 


±0,01 




110 


so/a 


0.350 


19.45 


2.70 


0.102 


18,36 


0,419 


20,12 


0,51 


0,31 


0.26 


0,23 


0.69 






±0.059 


±0.91 




±0.063 


±0.95 


±0.014 


±0.16 


±0,02 


±0,03 


±0.04 


±0.02 




129 


so/a 


0.489 


19.93 


2.70 


0.726 


20.95 


0.047 


16.68 


0.48 


0.65 


0.73 


0.82 


0.27 






±0.035 


±0.42 




±0.059 


±0.14 


±0.001 


±0.03 


±0,02 


±0,01 


±0.01 


±0,01 




135 


SO 


0.418 


20.45 


2.47 


0.490 


20,87 


0,039 


18,01 


0,76 


0,85 


0,88 


0,92 


0.43 






±0.046 


±0.69 




±0.048 


±0,20 


±0,009 


±0,37 


±0,05 


±0,03 


±0,02 


±0,01 




142 


so/a 


0.834 


21.40 


2.48 


0.318 


20,47 


1,094 


22,44 


0,84 


0,64 


0,39 


0,34 


0.36 






±0.150 


±0.85 




±0.127 


±0,55 


±0,072 


±0,25 


±0,01 


±0,01 


±0,02 


±0,02 




164 


so/a 


2.251 


22.62 


2.54 


2.541 


22,92 


0,061 


17,77 


0,46 


0,68 


0,91 


0,94 


0,12 






±0.188 


±0.45 




±0.207 


±0,13 


±0,007 


±0,19 


±0,06 


±0,03 


±0,00 


±0,01 




182 


SO 


0.355 


20.46 


2.40 


0.392 


20,77 


0,220 


21,91 


0,87 


0,85 


0.86 


0,90 


0,43 






±0.024 


±0.34 




±0,023 


±0,10 


±0,034 


±0,30 


±0,02 


±0,02 


±0,02 


±0,02 




209 


Sa 


0.625 


20.65 


2.59 


0.960 


22,06 


0,427 


20,93 


0,59 


0,50 


0.52 


0,64 


0.69 






±0.085 


±0.60 




±0,166 


±0,24 


±0,030 


±0,11 


±0,01 


±0,01 


±0,01 


±0,01 




211 


SO 


0.448 


20.12 


2.61 


0.458 


20.41 


0.432 


21.68 


0.85 


0.77 


0.74 


0.78 


0.61 






±0.051 


±0.46 




±0.051 


±0.15 


±0.052 


±0.18 


±0.01 


±0.02 


±0.03 


±0.02 




212 


E 


0.683 


20.78 


2.59 


0,781 


21,17 


0,052 


18,03 


0,70 


0,83 


0.89 


0,93 


0.09 






±0.036 


±0.32 




±0,034 


±0,09 


±0,009 


±0,27 


±0,04 


±0,02 


±0,01 


±0,01 




215 


SO 


0.457 


20.66 


2.45 


0,700 


21,77 


0,037 


16,95 


0,47 


0,64 


0.71 


0,81 


0.29 






±0.036 


±0.46 




±0,056 


±0,15 


±0,003 


±0,13 


±0,03 


±0,01 


±0,01 


±0,01 




233 


E/SO 


0.784 


20.40 


2.76 


0,952 


20,94 


0,113 


18,74 


0,64 


0,76 


0,86 


0,90 


0.35 






±0.086 


±0.76 




±0.104 


±0,23 


±0,014 


±0,17 


±0,08 


±0,05 


±0,03 


±0,01 




234 


Sb 


0.9% 


21.86 


2.42 


0.035 


17,80 


1,049 


22,03 


0,55 


0,22 


0,06 


0,05 


0.40 






±0.113 


±0.51 




±0.041 


±2.22 


±0.058 


±0.10 


±0,01 


±0,01 


±0.01 


±0,02 




236 


SO 


0.776 


21.10 


2.55 


0.944 


21.63 


0.055 


17.87 


0.59 


0.75 


0.86 


0.90 


0.29 






±0.047 


±0.37 




±0.055 


±0.11 


±0.004 


±0.12 


±0.04 


±0.02 


±0.01 


±0.01 




242 


E 


1.529 


21.92 


2.59 


1.825 


22.42 


0.263 


20.72 


0.69 


0.70 


0.87 


0.91 


0.10 






±0.171 


±0.76 




±0.214 


±0.22 


±0.015 


±0.24 


±0.07 


±0.06 


±0.02 


±0.02 




256 


E 


1.380 


20.88 


2.86 


1.551 


21.20 


0.090 


18.00 


0.63 


0.78 


0.91 


0.94 


0.14 






±0.073 


±0.33 




±0.071 


±0.09 


±0.014 


±0.25 


±0,05 


±0,03 


±0.01 


±0,01 




269 


HSO 


0.962 


20.39 


2.85 


1,045 


20,62 


0,093 


18,71 


0,78 


0,87 


0,94 


0,96 


0.35 






±0.058 


±0.39 




±0,056 


±0,11 


±0,019 


±0,28 


±0,05 


±0,03 


±0,01 


±0,01 




292 


so/a 


0.486 


20.32 


2.58 


0,297 


20,39 


0,56 1 


21,10 


0,64 


0,44 


0.36 


0,35 


0.26 






±0.137 


±1.35 




±0,218 


±0,93 


±0,022 


±0,26 


±0,02 


±0,03 


±0,04 


±0,01 




298 


SO 


0.554 


19.69 


2.83 


0,835 


20,84 


0,240 


19,53 


0,61 


0,60 


0,68 


0,78 


0.67 






±0.052 


±0.56 




±0.088 


±0,19 


±0,013 


±0,14 


±0,03 


±0,02 


±0,01 


±0,01 




300 


SO 


0.453 


19.84 


2.69 


0,584 


20,53 


0,092 


18,66 


0,66 


0,77 


0.82 


0,88 


0.41 






±0.028 


±0.43 




±0.039 


±0.13 


±0.002 


±0.10 


±0,03 


±0.02 


±0.01 


±0,01 




303 


E 


0.638 


20.59 


2.61 


0.774 


21.11 


0.091 


18.89 


0.67 


0.78 


0.86 


0.90 


0.15 






±0.058 


±0.63 




±0.068 


±0.19 


±0.012 


±0.16 


±0,06 


±0,04 


±0.03 


±0,01 




309 


HSO 


0.502 


20.02 


2.68 


0,540 


20,21 


0,060 


18,91 


0,87 


0,92 


0,94 


0,96 


0.14 






±0.023 


±0.31 




±0,020 


±0,08 


±0,016 


±0,40 


±0,04 


±0,02 


±0,02 


±0,01 




328 


Sa 


0.666 


20.72 


2,59 


0,862 


21,53 


0,392 


21,30 


0,76 


0,69 


0,71 


0,80 


0.47 






±0.082 


±0.58 




±0,108 


±0,19 


±0,046 


±0,20 


±0,02 


±0,02 


±0,02 


±0,01 




335 


so/a 


0.539 


20.86 


2.46 


0,541 


20,90 


0,059 


20,61 


0,94 


0,97 


0.98 


0,99 


0.37 






±0.040 


±0.51 




±0.033 


±0.14 


±0,053 


±1,28 


±0,06 


±0,04 


±0,02 


±0,01 




343 


SO 


0.404 


20.54 


2.43 


0,492 


21,10 


0,165 


20,93 


0,77 


0,80 


0,83 


0,88 


0.44 






±0.038 


±0.54 




±0.045 


±0.16 


±0.020 


±0.30 


±0,03 


±0,02 


±0.02 


±0,02 




353 


wso 


1.120 


21.05 


2.72 


1.205 


21.24 


0.134 


19.97 


0.81 


0,88 


0.94 


0.96 


0.02 






±0.051 


±0.32 




±0.050 


±0.09 


±0.021 


±0.22 


±0,04 


±0,03 


±0.01 


±0.01 




356 


Sa 


0.604 


20.20 


2.71 


0.392 


19.78 


0.899 


22.07 


0.90 


0.79 


0.65 


0.61 


0.46 






±0.044 


±0.31 




±0.035 


±0.11 


±0.039 


±0.17 


±0.00 


±0.00 


±0.01 


±0.02 




359 


SO 


0.574 


20.26 


2.67 


0,769 


21,05 


0,122 


19,05 


0,61 


0,71 


0.79 


0,86 


0.20 






±0.056 


±0.71 




±0,082 


±0,22 


±0,005 


±0,13 


±0,06 


±0,04 


±0,03 


±0,01 




360 


E 


0.341 


20.22 


2.45 


0,594 


21,66 


0,050 


17,62 


0,45 


0,61 


0,65 


0,77 


0,31 






±0.056 


±1.00 




±0,103 


±0,34 


±0,007 


±0,21 


±0,05 


±0,03 


±0,02 


±0,01 




366 


so/a 


0.630 


20.79 


2.55 


0,718 


21,16 


0,140 


20,46 


0,80 


0,85 


0,90 


0,93 


0.28 






±0.038 


±0.44 




±0.044 


±0,12 


±0,008 


±0,25 


±0,04 


±0,03 


±0,02 


±0,01 




368 


Sab 


0.542 


21.10 


2.39 


0.627 


22,39 


0,514 


21,56 


0,53 


0,39 


0,34 


0,41 


0.63 






±0.139 


±0.97 




±0.317 


±0.62 


±0.043 


±0.10 


±0,01 


±0,01 


±0.02 


±0,02 




369 


so/a 


1.036 


21.86 


2.44 


1.085 


21.98 


0.079 


20.28 


0.85 


0.92 


0.96 


0.98 


0.30 






±0.077 


±0.48 




±0.075 


±0.13 


±0.024 


±0.41 


±0.06 


±0.03 


±0.01 


±0.01 




371 


Sa 


0.956 


21.04 


2.66 


0.468 


20.64 


1.169 


21.92 


0.81 


0.63 


0.37 


0.34 


0.54 






±0.118 


±0.61 




±0.147 


±0.42 


±0.029 


±0.14 


±0.00 


±0.01 


±0.02 


±0.01 




372 


Sa 


1.435 


21.92 


2.56 


1.905 


22.72 


0.593 


21.99 


0.78 


0.71 


0.77 


0.84 


0.77 






±0.211 


±0.79 




±0.280 


±0.24 


±0.088 


±0.32 


±0,03 


±0,05 


±0.03 


±0,02 




375 


E 


4.979 


22.86 


2.82 


5,267 


23,01 


0,162 


19,22 


0,54 


0,69 


0,96 


0,97 


0.05 






±0.313 


±0.37 




±0,295 


±0,10 


±0,046 


±0,47 


±0,13 


±0,07 


±0,01 


±0,01 




381 


E/SO 


0.512 


20.34 


2.60 


0,550 


20,50 


0,088 


20,02 


0,88 


0,92 


0,94 


0,96 


0.37 






±0.063 


±0.89 




±0,060 


±0,25 


±0,041 


±0,64 


±0,10 


±0,07 


±0,05 


±0,02 




397 


so/a 


0.409 


20.50 


2.45 


0,111 


18,96 


0,540 


21,45 


0,68 


0,44 


0,35 


0,29 


0,57 






±0.035 


±0.45 




±0,026 


±0,36 


±0,014 


±0,10 


±0,01 


±0,01 


±0,01 


±0,01 




408 


SO 


0.382 


19.82 


2.63 


0,523 


20,66 


0,096 


18,86 


0,63 


0,74 


0.77 


0,85 


0.54 






±0.020 


±0.37 




±0.030 


±0.12 


±0.002 


±0.08 


±0,02 


±0,01 


±0.01 


±0.01 




409 


E 


0.498 


21.27 


2.30 


0.542 


21.49 


0.223 


22.60 


0.90 


0.90 


0.92 


0.94 


0.10 






±0.028 


±0.41 




±0.026 


±0.11 


±0.041 


±0.62 


±0,04 


±0,03 


±0.02 


±0,02 




410 


SO 


0.488 


20.67 


2.47 


0.575 


21,13 


0,157 


20,84 


0,80 


0,83 


0,87 


0,91 


0.43 






±0.032 


±0.46 




±0.037 


±0,13 


±0,013 


±0,29 


±0,04 


±0,03 


±0,02 


±0,02 




412 


E 


0.767 


21.39 


2.45 


1,091 


22,32 


0,056 


17,66 


0,42 


0,61 


0,76 


0,84 


0.32 






±0.072 


±0.56 




±0,096 


±0,17 


±0,007 


±0,22 


±0,06 


±0,03 


±0,01 


±0,01 




433 


SO/a 


1.021 


21.61 


2.51 


1,397 


22,46 


0,072 


17,95 


0,40 


0,59 


0.78 


0,85 


0.11 






±0.189 


±1.10 




±0,251 


±0,34 


±0,016 


±0,36 


±0,11 


±0,06 


±0,03 


±0,01 




440 


so/a 


1.108 


22.45 


2.29 


1,328 


22,98 


0,061 


18,80 


0,54 


0,72 


0.87 


0,91 


0.34 






























454 


so/a 


0.595 


20.42 


2.64 


0,990 


22,27 


0,461 


20,61 


0,47 


0,37 


0.37 


0,50 


0.58 






±0.068 


±0.50 




±0.182 


±0.26 


±0.020 


±0.07 


±0,00 


±0,00 


±0.00 


±0,01 




463 


SO 


0.639 


20.43 


2.66 


0.714 


20.74 


0.229 


21.18 


0.87 


0.87 


0.91 


0.94 


0.48 






±0.044 


±0.46 






±0.13 


±0.033 


±0.46 


±0.04 


±0.03 








465 


Sa 


0.609 


21.01 


2.47 


0.660 


21,22 


0,133 


21,10 


0,87 


0,90 


0.94 


0,96 


0.48 






±0.050 


±0.64 




±0.053 


±0,18 


±0,019 


±0,59 


±0,07 


±0,05 


±0,03 


±0,02 




481 


SO 


0.260 


20.25 


2.33 


0.157 


20,76 


0,280 


20,65 


0,32 


0,22 


0.22 


0,22 


0,74 






±0.036 


±0.65 




±0.090 


±0,70 


±0,004 


±0,12 


±0,03 


±0,04 


±0,04 


±0,02 




493 


E/SO 


0.193 


19.98 


2.28 


0,374 


21,66 


0,059 


18,74 


0,48 


0,61 


0,58 


0,73 


0.44 






±0.049 


±1.85 




±0,119 


±0,68 


±0,005 


±0,17 


±0,08 


±0,05 


±0,05 


±0,01 




523 


so/a 


1.044 


21.64 


2.51 


1,171 


21,96 


0,197 


21,08 


0,80 


0,83 


0.91 


0,94 


0.25 






±0.079 


±0.52 




±0,089 


±0,15 


±0.013 


±0.30 


±0,05 


±0,04 


±0.02 


±0,02 




531 


E 


1.549 


21.26 


2.80 


1.830 


21.71 


0.182 


19.30 


0.61 


0.69 


0.88 


0.92 


0.25 






±0.088 


±0.39 




±0.102 


±0.11 


±0.013 


±0.10 


±0.04 


±0.04 


±0.01 


±0.01 




534 


E 


0.620 


21.21 


2.41 


0.878 


22.17 


0.032 


16.68 


0.44 


0.63 


0.74 


0.83 


0.24 






±0.043 


±0.38 




±0.060 


±0.12 


±0.003 


±0.15 


±0.02 


±0.01 


±0.01 


±0.01 




536 


E 


1.266 


21.18 


2.73 


1.433 


21.52 


0.096 


18.57 


0.64 


0.78 


0.91 


0.94 


0.27 






±0.063 


±0.31 




±0.064 


±0.09 


±0.012 


±0.19 


±0.04 


±0.03 


±0.01 


±0.01 




549 


Sab 


0.866 


21.92 


2.34 


0.142 


19.90 


1.028 


22.50 


0.74 


0.47 


0.21 


0.17 


0.16 






±0.120 


±0.70 




±0.080 


±0.92 


±0.049 


±0.15 


±0.01 


±0.02 


±0.03 


±0.01 




Note. 


— Morphological types taken from Fabricatrt el at. (1999); r 


^slframc (fi-V)^ 


colors lakon from van Dokkum e 


al. (1998a); radii expressed in arcsec and surlact. 


brighmossos i 


1 mag per square 


arcsec. 



except for the fundamental plane parameters in which the radii have been transformed tokpc and the surface brightnesses toZ,Q/pc ; these structural parameters were derived by fitting the integrated profiles between 
r = o".l andr fa "^r-^ji- Apparent ellipticities derived from direct image fitting and are luminosity-weighted mean values for e. 
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Fig. 1. — Color images of the fundamental plane sample of CL1358+62 galaxies generated from the F606W and F814W HST imaging. The boxes are 15" on 
each side. The galaxies classified as E+A are ID # 209, 328, and 343. Galaxy # 234 is a star-forming emission line galaxy. Galaxy # 375 is the BCG. 
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Wavelength (A) 

Fig. 2. — The transmission curves for ttie WFPC2 filters F606W and F814W are superimposed witli tlie Johnson B, and Vj transmission curves redsliifted to 
z = 0.33. 
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Fig. 3. — The F814W .surface brightness profiles and integrated profiles from r = 1 pixel to ~2re. The thick solid lines show the fit of, and residuals from, the fit of 
the r'''*-law. The thin solid lines show the fit of and residuals from the bulge-plus-disk superposition, with the two components of the fit shown as dashed (r' ''-law 
bulge) and dotted (exponential disk) lines. The integrated curves were fit between r = i pixel and 2ri^2' Several galaxies are not well fit by a pure de Vaucouleurs 

profile, though as noted in the text, the effect of such deviations on the determinations of the fundamental plane parameter, re{I)g'^^, is quite small. The profiles and 
growth curves are plotted over arange of 0".l <r< 5", except for galaxies # 353, 371, and 375. 
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Fig. 3. (continued) ■ 
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Fig. 3. (continued) — 
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Fig. 3. (continued) — 
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Fig. 3. (continued) — 
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Fig. 4. — F814W residuals from the r'/^-law image modeling. Note the disky residuals in many galaxies. The boxes are 15" on each side. The galaxies classified 
as E+A are ID # 209, 328, and 343. Galaxy # 234 is a star-forming emission line galaxy. Galaxy # 375 is the BCG. 
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Fig. 5. — (a) The distribution of tlie ratio of half-Iiglit radii, derived from tlie best-fit r'/"-laws, to the half-light radii from the de Vaucouleurs (n = 4) fit. (b) The 
distribution of the ratio of half-Ught radii, derived from the bulge-plus-disk fits, to the half-light radii from the de Vaucouleurs fit. 
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Fig. 6. — (a) The ratio of half-light radii derived from all ;? ^4 to the effective radius from the de Vaucouleurs (n = 4) fit. (b) The ratio of the effective fundamental 
plane parameter from these n fits, with the effective fundamental plane parameter from the de Vaucouleurs fit. (c) The ratio of the mean effective fundamental plane 
parameter from these n fits, with the mean effective fundamental plane parameter from the de Vaucouleurs fit. The data for all 53 galaxies are shown, for every value 
of n fit. 
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Fig. 7. — (a-c) Bulge fraction, n, and ellipticity are plotted again,st the visual classifications from Fabricant et al. (1999). Note that, as expected, visually classified 
early-types have high bulge fractions and typically have large values of n. In (d) and (e) we plot the bulge fraction and n against ellipticity. In (f) we compare n with 
BF and note that those galaxies which are best-fit using large values of n tend to have high bulge fractions. 
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Fig. 8. — For every galaxy in the sample, we .show the ratios of {I)x / (/)4 vs. rxjr^, where X denotes either the (-'/"-law structural parameters (« = 1 ... 6) or the 
bulge-plus-disk structural parameters. Note that = rn,. For this plot, the structural parameters derived from every tested fitting range arc used, thus several points 
exist for each galaxy and profile shape. The soUd line indicates the curve of constant r(/) running parallel to the fundamental plane. The dashed line is the curve 
expected for a pure rV'*-law growth curve. 
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